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We study entropy generation in a one-dimensional (ID) model of bosons in an optical lattice expe- 
riencing two-particle losses. Such heating is a major impediment to observing exotic low temperature 
states, and "simulating" condensed matter systems. Developing intuition through numerical sim- 
ulations, we present a simple empirical model for the entropy produced in this ID setting. We 
also explore the time evolution of one and two particle correlation functions, showing that they are 
robust against two-particle loss. Because of this robustness, induced two-body losses can be used as 
a probe of short range magnetic correlations. 
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I. INTRODUCTION 



As cold gas experimentalists turn their attention to- 
ward more strongly correlated states (such as Mott in- 
sulators and fractional quantum Hall states) they find 
that equilibration times become large. This is problem- 
atic: inelastic processes limit the time over which one 
can conduct an experiment, and the experimentalists find 
themselves in a race. Can the system equilibrate before 
the entropy generated by the inelastic processes destroy 
the state of interest? A similar issue arises when the 
experimental protocol involves some sort of "adiabatic" 
change of parameters (such as ramping up the intensity 
of an optical lattice): can the adiabatic ramp be com- 
pleted before inelastic processes take over? Here we ex- 
plore a simple one-dimensional (ID) model where one 
can quantitatively study the entropy generated by two- 
body losses. Within this model we find that the entropy 
produced by each recombination event is of order the log- 
arithm of the number of atoms, highlighting the difficulty 
faced by future experiments. We hope that quantita- 
tive studies of such inelastic processes can help overcome 
them - though we see no simple solution at this point 
and time. On a more positive note, we find that within 
our model the inter-atomic correlations are largely time 
independent, even in the presence of drastic atom loss. 
Thus even with atom loss one can confidently measure 
the inter-atomic correlations of an initial state of interest. 
While one hopes that this result is generic, it is possible 
that the robustness of the correlations is an artifact of 
ID, where the dynamics are non-ergodic. To test this we 
have looked at small 2D clusters, finding nearly identical 
results. Extrapolating to more complex systems, we give 
an explicit example of how one could use losses to mea- 
sure nearest neighbor antiferromagnetic correlations in a 
two component Fermi gas. 

The problem of heating due to atom losses goes back 



to the first attempts to experimentally create degenerate 
Fermi gases 0, Q ■ Then it was pointed out that while 
atom losses are relatively benign in a Bose condensate, 
they have drastic consequences for a degenerate Fermi 
gas. Randomly removing particles from a Fermi sea gen- 
erates large amounts of entropy. 

Here we consider a version of a model introduced by 
Verstraete et al. Q. This and related models were ex- 
plored in a number of theoretical and experimental works 
mill. The original model consists of a gas of bosons 
moving in a ID lattice. Whenever two bosons are on 
the same site they recombine with rate Tq. The com- 
posite object which they form is then lost. This model 
could describe a gas of molecules (where a recombination 
mechanism always exists) or a gas of atoms (where light 
assisted collisions provide a recombination mechanism). 
Most experiments are engineered to minimize these two- 
body losses. They can itentionally be made stronger [l2[ 
and also occur in near-resonant optical traps, such as the 
blue-detuned lattices used by Schneider et al. [l3j . 

When To is sufficiently high, or in the presence of 
strong on-site interactions, one can integrate out the dou- 
bly occupied sites, producing a model with hard-core in- 
teractions and a nearest neighbor loss term. Mathemati- 
cally the time evolution of this dissipative system is then 
given by a Master equation 



dp 



(1) 



Qiajpa^al 



(riirijp + priiUj) 



where the conservative part of the dynamics are described 
by the Hamiltonian 



{i,3) 



a l a j + a ] a i 



(2) 
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In these equations, is the operator which annihilates 
an atom at site i, J is the tunneling matrix element, p is 
the density matrix, h is the reduced Planck's constant, 
(i, j) denotes nearest neighbor sites, 1/r is the time it 
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FIG. 1: (Color online) Left: Average particle number (TV) = ^2 n nTr p( n > as a function of time for an initial Mott insulator 
state on a L = 10 lattice. Solid line: numerical simulation; Dotted line: two-body decay law for an uncorrelated state 
N(t) = N(Q)/ [1 + 2Tt]. Middle: Same, but for a Tonks- Girardeau gas initial state (ground state of a hard core lattice gas 
with L = 10, iV = 6). Solid line: simulation. Green dashed curve: two-body decay law for an uncorrelated state. N(t) = 
N(0)/(1 + 2Tt), Dotted line: two-body decay low assuming time independent correlations N(t) = AT(0)/(1 + 2g (2) (O)n(O)R). 
Right: Average particle numbers in the different sectors (N^ n '(t)} = nTip^ n '(t) for the Mott insulator initial state. The sum 
of all curves at a certain time gives the blue in the leftmost figure. All times measured in units of the inverse hopping J . 



takes for atoms at nearest neighbor sites to recombine. 
These equations implicitly assume that one works in a 
Hilbert space where each site is occupied by only zero 
or one particle. The remarkable result found in previous 
work is that T scales inversely IV Thus when To — > oo, 
the dynamics become conservative, coinciding with those 
of a hard-core gas of particles, a lattice Tonks- Girardeau 
gas. 

We are interested in how entropy and correlations 
evolve with time in this model when r -C J, but T 0. 
For example, suppose one begins at time t — in an 
iV-particle Tonks state. As time evolves atoms are lost, 
until the system has N' < N particles. Does the sys- 
tem adiabatically evolve into a N' particle Tonks gas (or 
an ensemble of such gases with a different particle num- 
bers)? We find that this is not the case. As we describe 
below, we find that the N' particle system is better de- 
scribed by the initial iV-particle Tonks state with random 
atoms removed. This is a high entropy state. Despite 
its highly non-equilibrium character, it inherits the two- 
particle correlations of the initial A^-particle Tonks gas. 
These correlations, which are very different from what 
one expects for the AT'-particle Tonks gas, are directly 
measurable, and greatly impact the behavior of the sys- 
tem. 



II. NUMERICAL APPROACH 

To solve the Master equation in Eq. (fTJ), we consider 
a small chain of length L. We numerate all 2 L possible 
many-body states in which no more than one particle sits 
on each site. We explicitly write the density matrix in 
this basis as a 2 L x 2 L matrix, and express Eq. (fT]) as a 
coupled system of equations for the 2 2L matrix elements. 
Note the time evolution does not create coherences be- 
tween states of different particle number. This lack of 



coherence represents the fact that recombination events 
at different places or time are in principle distinguish- 
able. Consequently the density matrix is block diagonal, 
and for even number of particles N can be written 



P 



(o) 



P {2) © ■ 



P 



(N) 



(3) 



where p^ is the (JA -dimensional density matrix for the 
sector with n particles. We use a split-step method for 
our time evolution — alternating the exact Hamiltonian 
dynamics with the exact dissipative dynamics. We vary 
our time step to verify that our results are independent 
of the time step. We find that it is impractical to take 
L > 12, as the Hilbert space becomes too large. We 
typically quote results using L = 10orL = 12. 

Other approaches, such as the density matrix rcnor- 
malization group can be applied to this problem Q, al- 
lowing one to consider larger systems, but making some 
observables more difficult to calculate 



III. TIME EVOLUTION OF THE DENSITY 

As described by Garcia- Ripoll et al. [6J, the time evo- 
lution of the density in this model can be understood by 
a rather simple argument. One begins by noting that 
the rate of change of the number of atoms at a site only 
depends on the correlations between particles on nearby 
sites: 



d(nj) 
dt 



-Y [(nim+x) + (mru-i)] 



(4) 



Translational invariance implies that the two terms in 
brackets are equal to one-another, and one can write 



dn 
~d~t 



-2Tg^n 2 {t) 



(5) 
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FIG. 2: Time evolution of correlation functions starting from (left) the 10 particle Mott Insulator (L = 10, N = 10) or (right) 
the 6 particle Tonks-Girardeau state (L = 10, N = 6). Thick line: t = 0; Dashed line: t = 200J -1 ; Thin lines: intermediate 
times separated by 20 J -1 ; Dotted line: The single particle density matrix (a\ai+j) one would expect if each of the n-particle 
sectors were in their ground state at t = 200 J" 1 . The insets of the lower-left and lower-right figures show </ 2 ' as a function 
of density n = N/L together with the analytic formula for an infinite hardcore boson system in the ground state at the same 
density ffeq (n) = 1 — [sin(7rn)/(n7r)] 2 . 



where g^ = (n i n i+1 )/(n i ) 2 measures pair correlations 
and is related to the probability for finding two particles 
on neighboring sites. (Note: this differs from the def- 
inition of g^ in [13].) For uncorrected sites, such as 
one finds in the Mott insulator, one has g^ = 1. As 
we will directly illustrate in Sec. IIV1 when n < 1, the 
exact ground state of H has g^ < 1, and these equi- 
librium correlations are strongly number dependent. De- 
spite this behavior of the equilibrium correlations, the 
time evolution in Eq. ([1]) leaves the initial correlations 
nearly unchanged: <r 2 '(£) ~ 9^ 2 \t = 0). The degree to 
which this holds will be quantified below. If one treats 
the correlations as static one finds 

n[t) = l + 2trn(0) 5 ( 2 )(0)- ( } 

Variations on this equation have appeared in the litera- 
ture M, [l(| , which can be interpreted as different models 
for g®(t). 

Figure Q] (a, b) compares the time evolution found from 
the Master equation with that from Eq. ((6]). Panel (a) 
shows the evolution beginning from a Mott state (with 
g^(t = 0) = 1), while panel (b) shows the decay of a 
correlated state. Panel (c) shows the weight in each of 
the different particle number channels. 

Similar results for the decay of the Mott insulator were 
found in [Tfjj . In particular, Garcia- Ripoll et al. [3] used 
a Density Matrix Renormalization Group simulation to 



show that the observations in Fig. 1(a) are independent 
of system size. In the remainder of this paper we extend 
these results, directly analyzing the time dependence of 
correlations and entropy. 

IV. TIME EVOLUTION OF TWO-SITE 
OBESERVABLES 

We quantify the time evolution of the correlations by 
studying two objects: the single particle density matrix 
(ajaj+j) / (n) , and the density-density correlations func- 
tion (riirii+j) I (n) 2 . The equilibrium value of the lat- 
ter correlation function is efficiently calculated by per- 
forming a Jordan- Wigner transformation and mapping 
the hard-core Bose gas onto a gas of non-interacting 
Fermions. In particular, in the thermodynamic limit the 
equilibrium nearest neighbor density-density correlator 
at density n is 

roi , n ( sin 7m \ 2 , , 

g™(n) = l-( . 7 

4 \ nn J 

As already emphasized, we see large deviations from this 
equilibrium prediction. 

Figure [2] shows the time dependence of the correla- 
tion functions for the two initial conditions previously 
explored. The two figures on the left show the behavior 
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of the Mott state. Not only are the density correlations 
largely time independent, but so is the single-particle 
density matrix. One immediate implication is that the 
atom momentum distribution (and hence a time-of-flight 
image) is unchanged by the loss. The inset of the lower 
left figure compares the extracted from our simula- 
tions to Eq. ©. 

As shown by the two figures on the right, the dynamics 
from the Tonks state also leads to nearly time indepen- 
dent correlation functions. As time evolves there is a 
very slight drop in the nearest neighbor density corre- 
lations, and the single particle density matrix begins to 
fall off more rapidly with distance. This redistribution 
of the off-diagonal weight of the single particle density 
matrix corresponds to a shift of particles to larger mo- 
mentum. Interestingly, this is the opposite of what one 
would expect if one instead modeled the dynamics as just 
an adiabatic change in the number of particles. The in- 
set to the top-right figure shows the equilibrium single 
particle density matrix (corresponding to the number of 
particles are at time t = 200J -1 ). The slower slower 
spatial variation of the equilibrium (ajaj+j) / (n) , corre- 
sponds to a lower occupation of large k states. This is 
intuitively sensible, since when n = 1 one should have 
only the k = state occupied. 

In addition to being conceptually important, these cor- 
relation functions are directly observable. For example 
a time of flight measurement of the momentum distri- 
bution (as was for example done in [lij]) yields the 
Fourier transform of the single particle density matrix. 
The full density-density correlation function can be stud- 
ied within experiments with single-site resolution such as 
[l5j |. Alternatively noise correlation measurements [l6j 
or elastic light scattering pjj also probe this static struc- 
ture factor. 



V. ENTROPY 

We now proceed to calculate the entropy 

S(t) = -Trp(t)\np(t). (8) 

For ultra-cold atom experiments, where one is dealing 
with a small isolated system, the entropy is a more rele- 
vant than the temperature. This is especially true here, 
where the dynamics take one out of thermal equilibrium. 

It is convenient, as in Fig.[3l to parametrically plot S(t) 
as a function of N(t). For r< J the resulting curve is 
then independent of T. We find that the entropy is well 
approximated by the simple law 

S(J\0~lnW, (9) 

where Nq is the number of particles at t = 0. 

There is a particularly simple interpretation of this re- 
sult when one starts in the Mott insulating state. The 
entropy in Eq. Q is what one would find if one randomly 



punched holes in the Mott insulator. One would naively 
expect Eq. © to be an upper bound to the entropy, 
yet the numerical simulations find an entropy which is 
strictly above this curve. The extra entropy principally 
comes from the fact there is an indefinite number of par- 
ticles. When the Mott insulator is depleted to half-filling 
the entropy is maximal and one has an entropy per par- 
ticle of (for L, N — > oo) 

S/N ~ 21n(2) » 1.4. (10) 

A similar interpretation can be produced for the en- 
tropy of a depleted Tonks- Girardeau gas. There it is con- 
venient to map the system onto a gas of non-interacting 
fermions via a Jordan- Wigner transform. The fermions 
initially occupy No different momentum states. If one 
randomly removes fermions from momentum states one 
arrives at the entropy in Eq. This is similar to the 
arguments of Timmermans, where a trapped nonintcr- 
acting Fermi gas in three dimensions was considered Q. 
Again, the simulation finds an entropy slightly larger 
than this estimate. 

To understand the applicability of our results to higher 
dimensional systems, we have performed simulations of a 
two dimensional (2D) hardcore Bose gas on a rectangular 
square lattice geometry with a dimension of 4 x 3 sites and 
periodic boundary conditions. The result for entropy vs. 
particle number is also shown in Fig. [3] (a) and is almost 
identical to what we found in a ID system with the same 
number of lattice sites. 



VI. INDUCED LOSSES AS A PROBE OF 
LOCAL SPIN CORRELATIONS 

Our observation that initial correlations are preserved 
during time evolution makes two-body losses an ex- 
tremely powerful probe of cold atoms. Losses have long 
been recognized as a probe of correlations, but have al- 
ways been viewed as quite invasive [l(| El, E3-1HI • We 
find that two-body losses are a reliable way to learn about 
the initial correlations, even when a large fraction of the 
atoms are depleted. 

As a particular example, we now describe how pho- 
toassociation induced losses [IH, 0, [22| may be used to 
measure short-range magnetic correlations in two-species 
lattice bosons or fermions in the n = 1 Mott insulating 
phase 23]. This approach complements methods that 
detect long ran ge o rder, such as noise correlations and 
light scattering jl7l . |24| . Section IVI Al will discuss the 
fcrmionic case, while Sec. IVI Bl will deal with bosons. 

We wish to emphasize that this weak photassociation 
approach is very different from a sudden probe such as 
sweeping the system through a Feshbach resonance or 
photoassociating the system using a STIRAP protocol 
[13L [25j . In those latter approaches the signal size is 
limited by the instantaneous number of doubly occupied 
sites, which scales as J/U. For the weak probes used 
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FIG. 3: (Wide, Color Online) Left: Entropy S as a function of average particle number (N) during time evolution, starting 
from the (top, solid line) L = 12, N = 12 Mott insulator, (top, dashed-dotted line) 4 x 3, N = 12 (2D) Mott insulator, and 
(bottom, solid line) L = 12, iV = 6 Tonks-Girardeau initial states. Dashed line: analytic formula S ~ In (^j^) ■ Right: Entropy 
per particle as a function of time starting from the (solid line) L = 12, N = 12 Mott Insulator and (dashed line) L = 12, N = 6 
Tonks-Girardeau state with V = 0.01J. 



here, however, one can remove a large fraction of the 
atoms. 

The nearest neighbor s pin correlations are a smooth 
function of temperature |26j |. and represent an impor- 
tant precursor of the magnetic order which sets in on 
temperatures T ~ J s , where J s is the coupling constant 
in the effective spin model. Alternative approaches to 
measure these correlations use lattice modulation spec- 
troscopy [27| or manipulation of double well potentials 
combined with band mapping (28[. 



A. Two species fermions 



enough intensity that Jf ~ Jf- The subsequent density 
evolution will then be described by 

dn 
~dl 

where q is the number of nearest neighbors and the cor- 



(13) 



(2) 

relation function gl/ is given by 

gft = (rurij - 4S { ■ Sj}/ (riiY 



(14) 



The initial two-body loss coefficent is proportional to 

(2) 

gW(t = 0) = 1 — 4 (Si • Sj) and provides a direct measure 
of nearest neighbor spin correlations. 



We consider a two-species Fermi gas in a deep optical 
lattice that can be described by a Hubbard-model 



H 



-J E 



J7^n iit ni,4.(ll) 



We envision introducing a photoassociation laser which 
drives two atoms on the same site into a molecular state, 
which is lost from the system. The bare loss rate Tq, 
depends on the intensity of the photoassociation laser, 
as well as details of the atomic/molecular states. Inte- 
grating out the doubly-occupied sites produces a master 
equation similar to Eq. ([!]), which can formally be repre- 
sented as a a complex Hamiltonian 



H F = (j f + «4r) ( s > • s 



1 



(12) 



with J F = 4J 2 C/ / [Uo + (r /2) 2 )] and V = 
J 2 (r /2)/ [C/ 2 + (r /2) 2 ]. We envision letting the sys- 
tem equilibrate with the photoassociation lasers turned 
off (ro = 0). The lasers are then turned on at a low 



B. Two species bosons 

For an n = 1 Mott insulator of a two-species Bose 
gas one can similarly measure local spin correlation func- 
tions. Integrating out doubly occupied sites, the Hamil- 
tonian of such a system is formally (2t| 



Hb — ^ J z S*Sj — J± 



{SfS* + SfS]) 



(15) 



VniUj 



In the absence of losses, the coupling constants are re- 
lated to the hopping rate of the two species ( Jf , Jj_), and 
the on-site interactions between same (f/-M-, Uu) and dif- 
ferent (Ufi) species, 



J i = 



4J t J 



t J i. 



h = A 



U 



U- 



n 



Jy J+ 



2 


4 






4 " 






[u n 





V 



2U- 



1 





1 


1 




1.21% 





(16) 

(17) 
(18) 
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With a two species Bose gas one can se- 
lectively address three different photoassocia- 
tion transitions(tt^ m olecule, ti^molecule an d 
4-4-— > molecule) . This versatility may be used to measure 
a variety of nearest-neighbor spin-correlation functions. 
For example, driving a photoassociation resonance that 
converts an •f- and a 4--boson into a molecule can be 
formally described by substituting Uf± — > U^i — iT^/2. 
Specializing to the case — — J, corresponding to a 
typical optical lattice setup, one has 



H = [j*SfS? - Ji (SfSJ + SfS]) 



(19) 



+ o ( S i n i + niS i 



VriiTij 



-mnj - s*s* + S*S* + SfS^ 



where Jj_ - AJ 2 U n /[U^ + (r^/2) 2 ], J z = 
4J 2 E/ n /[L^ + (r^/2) 2 ] - 4J 2 (l/[/ tt + 1/U U ), V = 
J 2 /[U^ + (r n /2) 2 ] + J 2 (l/U n + 1/U U ) and T n = 
(r^/2) J 2 /[U^ + (r n /2) 2 ]. In this case two-body losses 



measure 



- (n^ - 4(S*S] 



5X QX 
i D j 



Sa&Mni) 2 (20) 



where i,j are nearest neighbors. Alternatively one 

(2) 

can photoassociate f"f— ^molecule to measure g^J = 



(rwrij + ASfS* + 2(S?nj + n l S'|))/(n l ) 2 . Photoassoci- 

(2) 

ating 4.4.— > molecule flips the sign of S z , giving g\/ — 
{mrij + 4S?S? - 2(Sfrij + n l 5|))/(n 4 ) 2 . Simultaneously 
photoassociating 4~t an d 4-4- a t the same rate allows mea- 
surement of (g<$ + gf^)/2 = (mrij + 45?5|). 

For a Mott insulating state of two species bosons, 
the expression in Ec|. ([20)) simplifies to g^(t — 0) = 
1 - 4(SfS? - S?S? - S?S?). Linear combinations of 
g$,g$\g$ probe (S.-S,) , (5?5?) and (5f5J + 5f^). 

Generically, for both bosons and fermions, the spin 
correlations at low temperatures tend to increase the loss 
rate. The ferromagnetic super-exchange in a Bose system 
encourages same-species atoms to sit next to each other, 
which due to Bose enhancement leads to an increased 
probability of doubly-occupying a site. The antiferro- 
magnetic super-exchange in a Fermi system enhances the 
probability of up-spins lying beside down-spins, increas- 
ing the chance that they will end up on the same site. 
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